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ABSTRACT 

In the mean-field theory of magnetic fields, turbulent transport, i.e. the turbulent electromotive force, is de- 
scribed by a combination of the a effect and turbulent magnetic diffusion, which are usually assumed to be 
proportional respectively to the mean field and its spatial derivatives. For a passive scalar there is just turbulent 
diffusion, where the mean flux of concentration depends on the gradient of the mean concentration. However, 
these proportionalities are approximations that are valid only if the mean field or the mean concentration vary 
slowly in time. Examples are presented where turbulent transport possesses memory, i.e. where it depends 
crucially on the past history of the mean field. Such effects are captured by replacing turbulent transport co- 
efficients with time integral kernels, resulting in transport coefficients that depend effectively on the frequency 
or the growth rate of the mean field itself. In this paper we perform numerical experiments to find the charac- 
teristic timescale (or memory length) of this effect as well as simple analytical models of the integral kernels 
in the case of passive scalar concentrations and kinematic dynamos. The integral kernels can then be used to 
find self-consistent growth or decay rates of the mean fields. In mean-field dynamos the growth rates and cycle 
periods based on steady state values of a effect and turbulent diffusivity can be quite different from the actual 
values. 

Subject headings: MHD - turbulence 



1. INTRODUCTION 

A simple form of turbulent transport is the mixing of a 
passive scalar associated with the mutual exchange of fluid 
parcels. This process is similar to non-turbulent mixing that 
occurs just because of thermal fluctuation or Brownian mo- 
tion, often referred to as molecular diffusion. The latter pro- 
cess is described by a diffusion equation with a diffusion term 
of the form kV 2 C, where k is the molecular diffusion coef- 
ficient and C is the concentration. Turbulent diffusion, on 
the other hand, applies to a suitably averaged mean concen- 
tration, C, and is normally described by a diffusion term of 
the form Ar t V 2 C, where /c t is a turbulent diffusivity. The ratio 
KxIk scales like the Reynolds number (or, more precisely, the 
Peclet number) and can become very large under many astro- 
physical conditions (stars, accretion discs, galaxies). 

Problems connected with this simple prescription occur 
when the mean concentration shows variations on timescales 
shorter than or comparable to the correlation time of the tur- 
bulence. In practice this means that a sinusoidal profile of C 
with wavenumber k would decay at a rate K t k 2 where K t is no 
longer constant, but it depends itself on the actual decay rate. 

The fact that problems occur when the mean concentration 
changes on short timescales should not be surprising. In- 
deed, in the textbooks of Moffatt (1978) and Krause & Radler 
(1980) it is shown that a proper description of turbulent trans- 
port involves a convolution of an integral kernel with the mean 
concentration over past times. This is why one talks about 
memory effects: the turbulent diffusion is not just an instanta- 
neous property of the turbulence, but depends on its full time 
history (Hori & Yoshida 2008). Dealing with a convolution 
over past times is an unpleasant complication, so its effects 
are often neglected. However, there can be circumstances of 
astrophysical relevance where this is no longer permissible. 
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Such a circumstance is the damping of solar p-mode oscil- 
lations through turbulent motions in the surface layers (Stix 
et al. 1993). Here the timescales of /5-modes and convection 
are comparable, so memory effects must be important. Stix et 
al. (1993) find that the turbulent diffusion is reduced by a fac- 
tor exp(-a»oscT), where <x> osc is the oscillation frequency and 
t is the correlation time of the turbulence. Memory effects 
have also been invoked in connection with propagating front 
solutions in the galactic dynamo (Fedotov et al. 2002, 2003), 
and variations of the solar cycle (Otmianowska-Mazur et al. 
1997), although there the timescales are more disparate. 

A practical way of dealing with memory effects has been 
proposed by Blackman & Field (2002, 2003), who derived 
an evolution equation for the turbulent flux of concentration 
based on a simple closure prescription known as the r approx- 
imation. One of the main beauties of this approach is that the 
usual diffusion equation, which is of parabolic nature, is now 
replaced by a damped wave equation, which is of hyperbolic 
nature. This implies that signal propagation is no longer in- 
finitely fast, but its speed is limited to the rms velocity of the 
turbulence. The principal validity of this approach has been 
demonstrated using turbulence simulations of passive scalar 
diffusion (Brandenburg et al. 2004). One of the goals of the 
present paper is to provide a more direct means of determin- 
ing memory effects of turbulent transport that can also be ap- 
plied to more complicated cases of vector fields such as the 
magnetic field. 

A promising method for calculating turbulent transport co- 
efficients for the magnetic field is the test-field method. In 
this approach one calculates evolution equations for the small- 
scale field that results from a given set of different test fields. 
In this way one can calculate the full tensorial nature of the 
turbulent diffusion tensor, as well as the a tensor that can be 
relevant for amplifying the magnetic field if the turbulence 
lacks mirror symmetry, for example in the presence of helic- 
ity. These test fields have a given length scale characterized 
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by some wavenumber. By varying this wavenumber it has 
been possible to determine the scale dependence of the mean 
fields that are being diffused and/or amplified (Brandenburg et 
al. 2008a). Using a Fourier transformation over all wavenum- 
bers, it is possible to determine the spatial properties of the 
integral kernels that are used in the convolution with the mean 
field over all other points in space. It is customary to approx- 
imate the kernels by 6 functions, in which case the convo- 
lutions become multiplications. In the test-field method, the 
corresponding coefficients are obtained as the limit of vanish- 
ing wavenumber. However, in order to make statements for 
finite domains of length L, the wavenumber k — k\ = 2n/L is 
most relevant. Unless stated otherwise, we focus therefore on 
results for k = k\. 

In an analogous fashion, we can make the test fields time- 
dependent and compute in this way the temporal properties of 
the integral kernels. By imposing sinusoidal variations of the 
test fields over a range of different frequencies we calculate 
the integral kernels first in Fourier space, because there the 
convolution corresponds just to a multiplication. The integral 
kernel in real space is then obtained by Fourier transforma- 
tion. Another possibility is to apply an exponentially growing 
or decaying time variation. In a sense this comes closest to the 
application of calculating modifications of growth rates due 
to finite memory effects. The integral kernel can then be cal- 
culated by inverse Laplace transformation, but this approach 
involves integration along the imaginary axis and is therefore 
only feasible if the data can be fitted reliably to an analytic 
function. We note that it is in principle also possible to de- 
termine integral kernels directly by applying a 5 function-like 
variation to the mean concentration gradient or the mean field, 
but the disadvantage here is that it is then not so easy to im- 
prove the statistics by time averaging. Nevertheless, such a 6 
function-like perturbation provides an additional verification 
and is certainly a useful thought experiment. 

The temporal properties of integral kernels in turbulent 
transport may be particularly important in dynamo theory 
where simulations and theory are now sufficiently accurate to 
show finite memory effects under controlled conditions. As 
a side effect, growth rates based on a dispersion relation with 
constant a effect and turbulent magnetic diffusivity may be- 
come inaccurate. It is quite plausible that under more compli- 
cated circumstances finite memory effects will be even more 
important. However, without proper knowledge of what to ex- 
pect, this would only remain speculation. A goal of this paper 
is therefore to clarify finite memory effects in simulations of 
forced helical turbulence in a periodic domain. We consider 
here only the kinematic case, i.e. the velocity is unaffected by 
the magnetic field. 

In Section|2] we will motivate our work by considering two 
approaches to calculating the growth rate of the Roberts flow 
dynamo. In Section [3] we define our formalism, most im- 
portantly the time response kernels that describe "memory" 
effects. We will treat both the turbulent transport of mag- 
netic fields and the conceptually simpler transport of passive 
scalars. In Section|4]we give a brief theoretical overview be- 
fore discussing our numerical simulations and results in Sec- 
tions [5] and [6] We discuss those results in Section [7] and con- 
clude in Section[8] 

2. BACKGROUND: MISMATCH IN GROWTH RATES 

A direct approach to determining the growth rate of a dy- 
namo is to solve the induction equation for the magnetic field 



B numerically: 

dB , 

— = V x (U x B) + r]V 2 B. (1) 

at 

Here U is the velocity and 77 is the microscopic magnetic diffu- 
sivity. We are interested in dynamos that produce mean fields, 
B, denoted here by an overbar. In the following we take this 
to be an xy average. We calculate then the growth rate of the 
mean field as _ 

/Igrowth = din fi rms /df. (2) 

This can now be compared with the corresponding result from 
mean-field theory, where one considers the averaged induc- 
tion equation, 

dB - - - 

— = Vx({7xfi + £) + nV 2 B, (3) 
at 

with _ 

£, = uxb (4) 

being the turbulent electromotive force and u — U - U and 
b = B -B are the fluctuations. Symmetry considerations con- 
strain the form of £, and in the case of homogeneous isotropic 
turbulence with helicity, the expression for £ is found to be 

8-aB- r/t/JoJ, (5) 

where a describes the a effect, j] t is the turbulent magnetic 
diffusivity, J = V x B/fio is the mean current density, is 
the vacuum permeability, and higher order terms have been 
omitted. Such a model is generally referred to as an a 2 dy- 
namo. For references see Moffatt (1978) and Krause & Radler 
(1980). 

A new and accurate method for determining a and rj t is the 
test-field method of Schrinner et al. (2005, 2007) that will be 
described below. The details of this method are not essential 
at this point, except that we do emphasize that for our values 
of the magnetic Reynolds number Re^ the wavenumber of 
the test field is chosen to be that of the box, which is also the 
smallest wavenumber that fits into the domain. 

For isotropic turbulence in a periodic domain the magnetic 
field can develop long wavelength variations in any of the 
three coordinate directions (Brandenburg 2001). We assume 
this to be the z direction and use averages over the x and 
y directions. Solutions of a homogeneous a 2 dynamo obey 
Vxfi = k : B = fioJ and are proportional to exp(ik,z + At) with 
the dispersion relation 

A = ak z + rj^k 2 , (6) 

where k, is the wavenumber in the z direction. Both a and rj t 
are taken as constant in space owing to the assumed statistical 
homogeneity of the turbulence. For flows with positive kinetic 
helicity, a is expected to be negative, so growing solutions 
correspond to negative values of k z . 

We refer to the value of A obtained from the dispersion re- 
lation (O as /tdisp. This is the second approach to determining 
the growth rate of the dynamo. It has the disadvantage of be- 
ing indirect, but the advantage of aiding comprehension of the 
dynamo mechanism itself. If the theory behind this second ap- 
proach is correct, then the results should match, so comparing 
the growth rates allows one to test the validity of Equation 

In order to motivate the purpose of this paper, let us now 
compare in Figure [T|/lg rowtri with A& sp for the simpler case of 
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Fig. 1. — ReM dependence of the growth rate for the Roberts flow as ob- 
tained from a direct calculation (/lgrowth) compared with the result of the dis- 
persion relation, Aa sv = ak z — (ij + il t )k}, using a cubic domain of size L? , 

where k\ = 2nfL and kf = y[2k\ . For this range of Re^, , the most unstable 
mode is the largest one that fits in the box (k, = k\ ). The two horizontal lines 
in gray mark the values of /ig low th ar| d ^disp at Re*/ = 30, denoted by (i) and 
(ii), respectively. 

a steady periodic helical flow instead of turbulence. We use 
here the Roberts flow, whose details will be discussed later. 
The two estimates for A do indeed agree when A = 0, at the 
critical magnetic Reynolds number for the onset of dynamo 
action Recent - 5.52. For larger values of Re^, there is a 
discrepancy that can become rather dramatic for Re^ > 20. 

One of the motivations for our work then is the fact that 
the growth rate estimated from Equation ©, where a and rj t 
are obtained from the test-field method, becomes increasingly 
inaccurate for large growth rates, implying that Equation © 
is inadequate to describe growing dynamos. We emphasize 
that this discrepancy vanishes not only in the marginal case, 
but also for the nonlinearly saturated dynamo. This is why in 
Brandenburg et al. (2008b) the quenched values of a(B) and 
r]i(B) were found to obey Equation (O with A — 0. 

Even though the Roberts flow has been studied extensively 
over the years (see, e.g., Roberts 1972, Soward 1987, Plunian 
et al. 1999, Plunian & Radler 2002a,b, Courvoisier 2008), and 
especially so in connection with the Karlsruhe dynamo exper- 
iment (cf. Stieglitz & Miiller 2001, Radler et al. 2002), a dis- 
crepancy between theoretically expected growth rates based 
on mean-field theory and the actual ones has never been re- 
ported. For example in Plunian & Radler (2002a), the ac- 
tual growth rates have been determined directly without in- 
voking mean-field theory, and in Radler et al. (2002) only the 
marginal case has been compared with observations. How- 
ever, in the marginal case the discrepancy disappears. In Plu- 
nian & Radler (2002b), on the other hand, the values of a and 
T] t have again been determined self-consistently for cases dif- 
ferent from the marginal one. Thus, the mean field is then of 
course no longer steady, and so their values of a and j] t apply 
only to this particular time dependence, but not to a fictitious 
steady case, for example. We say here "fictitious", because 
for given values of Re^ and k z , there is normally only one 
relevant solution, namely the one with the largest value of A. 
However, for a predictive theory one should know a and rj t 
before having solved the problem, i.e. before knowing A. In 
the following we explain how the fictitious steady case can 
actually be realized in a simulation. 

In order to clarify the point that, for given values of Re^, 
a and rj t depend also on the resulting growth rate, let us now 



Fig. 2. — Dependence of A on A for Rsm = 30. The values of A = /ldi sp 
and /lgrowth °f Figure [T] are indicated by a vertical dashed and dotted lines, 
denoted by (ii) and (i), respectively. Note that A = (dash-dotted line) for 
A = 

/Idisp* where /ldisp — a k z — (tj + i]t)k z with a and t]i obtained using the 
test-field method for steady fields. The linear interpolation between the points 
(A, A) = (0, igrowth) and (i gl 

owth> 0) is indicated by a triple-dot-dash line. 

consider a modified induction equation with an artificial "fric- 
tion" term, 

^=Vx([/xB) + ifJ 2 B - KB, (7) 
of 

where A is a new control parameter and B is the xy-averaged 
field. Note that the evolution of the departure from this xy- 
averaged field, b = B - B, is unaffected by this manipulation, 
so £ = m x ft is exactly the same as before. The solutions 
of B have still an exponential time dependence, and standard 
mean-field theory gives for the growth rate A of the mean field 

A = ak z -{r] + j] t )k 2 z - A. (8) 

So, as the value of A is increased (for given values of Re^ and 
k z ), the growth rate A decreases. [The tilde has been added 
to distinguish A from that used in Equation ([6j.] There is a 
critical value A* for which 2 = 0. This value is determined by 

A. = ak, - (rj + T] t )k 2 . (9) 

Given that in this case the mean field is steady, we now expect 
Equation (O to be accurate. To verify this we solve Equa- 
tion (0 numerically and determine the growth rate A. The re- 
sult is shown in Figure|2]where we plot A vs. A for Re^ = 30. 
For A = we find A = /l gr0 wth- More importantly, it turns out 
that A = at a value A = A* = /Idisp, indicated by (ii), that is 
given by Equation (O with the same values of a and rj t that led 
earlier to the discrepancy in Figure Q] Most crucially, the nu- 
merically determined growth rate A in Figure|2]deviates from 
a linear interpolation between the points (A, A) = (0, /lgrowth) 
and (/lgrowth, 0). This suggests again that the assumption of the 
a and rj t in Equation ([8]) being independent of A is invalid. 

We note that for larger values of Re^ (e.g. for Re^ = 50), 
Equation (Q permits additional solutions with insignificant B 
that cannot be damped by the AB term. However, as a proof 
of concept, it was only essential that Re M was big enough so 
that there is a clear difference between A gmw ± and /Idisp- 

The results presented above show that a naive application 
of the dispersion relation to cases where A + is not possible 
and gives results that disagree with the direct simulation. This 
is because the values of a and r\ x apply only to the steady case, 
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as demonstrated by considering the associated steady problem 
of Equation 0, where A = A, is predicted from Equation (0 
using the a and rj t values obtained from the test-field method. 

Recently, Hori & Yoshida (2008) noted that, in the Roberts 
flow, memory effects can be responsible for an enhancement 
of the growth rate. The reason why Plunian & Radler (2002b) 
found the correct growth rates from Equation (0 even when 
A + is that their values of a and t]^ were automatically 
"tuned" to the resulting growth rate. Their values do therefore 
not apply to the steady case, which can be verified by consid- 
ering the mean-field problem associated with Equation 0. 

To understand the reason for the discrepancy between ac- 
tual growth rates and those obtained from the standard (time- 
independent) test-field method, it is important to recall that 
a multiplicative relation in Equation 01 is only an approxi- 
mation and that it should instead be a convolution in space 
and time (Moffatt 1978; Krause & Radler 1980). Alterna- 
tively, a Taylor series expansion of £ in space and time can 
be employed. Already in the simple case of the Roberts flow 
Equation (0 cannot be justified when the mean field changes 
sufficiently rapidly in time. In this paper, we show that in such 
cases "memory" effects of the turbulent transport coefficients 
cannot be ignored. This implies that the electromotive force 
at a given time depends not only on the mean fields at that 
specific time, but also on the mean fields at all prior times. 
In practice, this means that the turbulent transport coefficients 
depend themselves on the resulting growth rate and/or fre- 
quency of the mean fields. 

3. FORMALISM 

Quite generally, we are interested in expressing quadratic 
correlations of fluctuating quantities in terms of mean fields. 
Examples include the mean turbulent concentration flux and 
the mean turbulent electromotive force, 

T = uc, and & = uxb, (10) 

respectively. Here, c — C - C is the fluctuation of the concen- 
tration density. The number of preferred directions available 

to mean quantities such as T and £ are limited, and so the 
aim is to relate them respectively to the gradient of the mean 

concentration, G = VC, and to a linear combination of the 
mean magnetic field B and its curl, V x B — (XoJ. However, 
instead of multiplicative (instantaneous) relations of the form 

T = -K t G, T, = a~B-r lx nJ, (11) 

we now adopt such relations in their more general forms in- 
volving a convolution in time, i.e. 

T(t) = - f k t (t-t')G(t')dt', (12) 

and 

6(0 = f &(t-t')B(t')dt'- f Ut-f)fio7(t')dt', (13) 

\J — CO vV —Of! 

where quantities with a hat denote integral kernels, so k t (t) 
is an integral kernel describing turbulent passive scalar diffu- 
sion, a(t) describes the a effect, and fj t {t) the turbulent mag- 
netic diffusion. This approach is the most general search for 
memory effects, and we adopt it to find out how to modify 
Equation (0 to model more accurately growing dynamos. 

We recall that in general our averages (being two dimen- 
sional over the xy plane) are also functions of z, but the z 



dependence has here been suppressed in favor of a more com- 
pact notation. In general, Equations (V2\ and dT3l > should also 
include a convolution over z- This property has recently been 
studied in Brandenburg et al. (2008a), but the spatial aspects 
of the convolution will here be ignored by considering mag- 
netic fields that have only a single wavenumber k z , which cor- 
responds to the smallest wavenumber k\ = 2n/L that fits into 
the domain of size L? . 

3.1. Standard test-field methods 

In this section we reiterate the essence of the standard 
test-field methods for calculating a, rj t , and k u where mem- 
ory effects are ignored. As noted above, mean-field theory 
treats turbulent transport through the correlations of fluctuat- 
ing quantities as in Equation (fTTb . If the transported quan- 
tity does not itself affect the dynamics of the system, as in 
the cases of passive scalars or kinematic dynamos (where the 
magnetic field is too weak to affect the momentum equation), 
then the transport coefficients are functions of the velocity 
fields alone. 

This independence of the transport coefficients of the mean 
field implies that the transport coefficients will be found also 
in systems where a mean field is externally imposed and does 
not obey any evolution equation. Such a field is called a test 
field. A set of different test fields is needed to_ determine si- 
multaneously the prefactors a and rj t of B and J, respectively. 
In the test-field method of Schrinner et al. (2005, 2007), one 
subtracts the mean-field equation from the full induction 
equation (0 to obtain an evolution equation for the fluctuating 
field b, 

db - - , 

— = Vx(Uxb + uxB + uxb-uxb) + rf! b. (14) 
dt 

This equation is then applied separately to each of the fields 
B pq , where p — 1 or 2 and q - cor s label different test fields. 
Brandenburg et al. (2008a,b) use the four test fields 

B lck = B (cosfe,0,0), B lsk = B (sinfe,0,0), (15) 

B 2ck = B (0, cos kz, 0), B 2sk = B (0, sin kz, 0), (16) 

where the third superscript k has been added to denote the 
wavenumber, and Bo is a normalization factor. The response 
to each test field, b pqk , is found by solving Equation (TBI . In 

this way, one finds &" lk - ux b pqk and obtains 4x2 equa- 
tions, 

■ppqk -jP4 k 

& \ = a U B j -ntjfJeJj > (17) 

for the 4 + 4 unknowns, a,j and 77,7, for i = 1,2 and j = 1,2. 
These eight unknowns are obtained as 

( a u \_ J cos kz sinfe \/£f*\ fl „ 
\mflk}~ l-sinfe aMltz)\e!• k )^ { ' 

where the rank-3 tensor 77,^3 is related to the rank-2 tensor in 
Equation ( fTTb via 77,7 = rjikiejis- Note that the result is inde- 
pendent of the value of Bq. For stationary isotropic homoge- 
neous turbulence we have constant values of a\ \ = an = a 
and 7711 = 7722 = 77 t , except for statistical fluctuations resulting 
from finite computational volumes. 

The test-field method for a passive scalar works analo- 
gously (Brandenburg et al. 2009). The concentration per unit 
volume C obeys the equation 

— = -V • (UC) + kV 2 C, (19) 
dt 
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and the evolution of the mean concentration C is obtained by 
averaging Equation $1% , which yields 



— = -V • (UC + T) + kV 2 C. 
at 



(20) 



The test scalar equation is obtained by subtracting Equa- 
tion d20b from Equation ( fT9] l, which yields 

S = - V • (Uc + uC + uc - uc) + kV 2 c. (21) 
at 

In order to obtain K t , one uses the test scalars 



(? = Co cos kz, C* = Co sin kz, 



(22) 



where q = c or s denotes the spatial dependence of the 
test scalar and, again, an additional superscript k denotes the 
wavenumber, while Co is a normalization factor. For each test 
scalar, we obtain a separate evolution equation for c qk . In this 

way, we calculate the fluxes, T qk - uc qk , and compute the 
three components of : 

KB =-{- sin kz¥f + cos kzWf) /kC , (23) 

for i = 1, 3, where ( ) z denotes a z average. Again, the val- 
ues of Ki3 are independent of the normalization constant Co- 
For stationary isotropic homogeneous turbulence we have, ex- 
cept for statistical fluctuations, constant Ky = K t Sij. 

By applying the test-field and test-scalar methods to a range 
of different wavenumbers k, it was possible to assemble two 
full integral kernels in space (Brandenburg et al. 2008b, 2009) 
and hence to take the effects of finite scale separation into ac- 
count. In the following, we proceed analogously by applying 
the test-field and test-scalar methods to a range of different 
frequencies to assemble two full integral kernels in time and 
hence to take memory effects into account. 

3.2. Determination of the kernels 

As is common in linear response theory, all integral ker- 
nels vanish for t < 0. Therefore the integrations in Equa- 
tions (fT2l and ( TT3l extend effectively only to t' = t. In order 
to determine these kernels numerically, we can either calcu- 
late them directly by imposing 6 function-like variations of 
the test fields, or we can use the fact that a convolution corre- 
sponds to a multiplication in spectral space, i.e. 



where 



T(oS) = -k t (to)G(a>), 



k t (oj) = dte i0J 'k t (t) 



(24) 
(25) 



is the Fourier transform of k t (t). 

A multiplicative relation between T and G applies also to 
the Laplace transform of these functions with 



where 



no = -* t (i)G(i), 



i<t(s) = I dte k t (t) 
Jo 



(26) 
(27) 



is now the Laplace transformation of k t (t). 

We introduce an additional superscript lo for the cases 
where the test fields or concentrations have cos tot time depen- 
dence and superscript s for the cases where the test fields or 



concentrations have exp st time dependence. The superscripts 
or the explicit time dependence are sometimes suppressed. In 
most of the cases we use test fields with a sinusoidal spatial 
dependence with wavenumber k — k\. However, it is some- 
times useful to vary also the value of k. In these cases, we also 
add the superscript k. 

The multiplicative relations above imply that for an oscil- 
latory perturbation with a single frequency there is a multi- 
plicative relation between G qkw (f) and < F qka> (t), where the first 
superscript denotes the frequency; see Appendix [Al In gen- 
eral, K t is a tensor, but in the following we restrict ourselves 
to determining only one of its components, namely the one 
relating the z components of 'Fit) and G(t) to each other. We 

therefore assume G(z,t) = (0, 0, G), where G = dC/dz. The 
different test scalars C qku> are denoted by superscripts c and 
s for spatial dependences proportional to cos kz and sin kz, so 
we have 

C = Co cos kz cos cot, C = Co sin kz cos cut, (28) 
for oscillatory test fields, and 

C = Co cos kz exp st, C — Co sin kz exp st, (29) 

for exponentially growing or decaying test fields. For each 
value of u> we determine the resulting z component of the flux, 

faj(t)- As shown in Equation ( lA4b of Appendix [Al we can 
calculate the response kernel as 



k t (k,oS) = -2G l (e la "T (0 



—koj 



(30) 



where the subscript t behind an angular bracket denotes a 
time average. Note that ki(k, a>) is complex such that its real 
part is symmetric about to — 0, while the imaginary part is 
antisymmetric. In other words, it obeys the Kramers rela- 
tion, k t (k, —cS) = k\(k, u>), where the asterisk denotes complex 
conjugation; see, e.g., Moffatt (1978) and Krause & Radler 
(1980). In our case, in addition, k t is symmetric in k. 

Analogous relations apply to !F(s) and G(s). In this case, 
Equation ( 1301 ) is modified to 



K t (k, s) 



-A-.v 



e~ a r it) 



(31) 



As discussed in Section [3TT1 our test fields allow us to pick 
out each tensor component of a t j and rjij separately. We there- 
fore define time-dependent test fields 

gpqkm = gpqk CQS ^ an£ j gpqks = gpqk ^ ^ ^ 

where the time-independent test fields B pqk were defined in 
Equation ( TTTb . Owing to variations of the form sinfe and 
cos kz one multiplies with the inverse of a rotation matrix, 



( &ij(k,s) 
\ fjij3(k, s)k 



)-H 



cos kz 
- sin kz 



sinfe 
cos kz 



(33) 



where the matrix above results from the choice of the sinu- 
soidal test fields; see Sur et al. (2008) for details. An analo- 
gous equation applies also to the case of oscillatory test fields 
where s is replaced by -\co, so we write 

/ &ij(k,aj) \_l m>t( cos kz sin kz \ I \\ r , 
\fj ij3 (k,cj)k)~\ e v -sinfe cosfe^£2^j/ ( ' '- 4i 



keeping in mind that a tilde has been used to indicate both 
Fourier and Laplace transformation. 
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4. PRELIMINARY CONSIDERATIONS 

Before entering the numerical determination of the integral 
kernels let us consider a current approach that captures mem- 
ory effects, as well as its simplest extension. This will later 
serve us with a useful fit formula for the more complicated 
cases. 

4.1. Expectations from the r approximation 

We use the term t approximation here in the form intro- 
duced by Blackman & Field (2002, 2003, 2004). The essence 
of the r approximation is to write down evolution equations 
for second order correlations such as uc and u xb. This 
results in triple correlations that are not omitted, as in the 
first-order smoothing approximation (FOSA), but are instead 
approximated by a closure hypothesis. In the t approxima- 
tion, one replaces the triple correlations by quadratic correla- 
tions divided by a relaxation time t (Vainshtein & Kitchatinov 
1983; Kleeorin et al. 1996). This timescale is expected to be 
comparable to the turnover time of the turbulence. 

Blackman & Field (2002, 2003, 2004) were the first to re- 
tain the time derivative in the evolution equations for uc and 
u x b. This means that the Fickian diffusion approximation 

of Equation ( fTTI ). i.e. T = -£ t ()G, with £ t o = tu 2 for one- 
dimensional diffusion in the z direction, is generalized to 



(35) 



This implies that the Fourier-transformed integral kernel is 



km 



Tilt 



1 — icot 1 - ilot 



(36) 



(Any k dependence is here ignored.) In real space, this ex- 
pression for K t (co) corresponds to the integral kernel 



C — 

J-co 2n 



Tilt 



1 — iu>T 



u 2 0(f) e 



-t/T 



(37) 



where the integral has been solved as a contour integral 
around the pole at to — -i/r, and © is the Heaviside step 
function with 0(f) = 1 for f > and otherwise. In the limit 
r — > 0, the exponential function reduces to T(5(f), so 



k t (t) -> tu 2 6(t) = Kto6(i) (for r -> 0), 



(38) 



and one recovers the usual prediction in which turbulent dif- 
fusion can be treated as a multiplicative enhanced diffusion 
coefficient. 

Similar considerations also apply to the case with magnetic 

fields, where £ is essentially being replaced by (1 + t8,)£>. For 
exponentially growing solutions, one would therefore expect 
that the actual growth rate A is reduced by a factor (1 + /It) -1 . 
However, this expectation may be too naive and will need to 
be reconsidered in this work. 

A useful diagnostic for the applicability of Equation d36l l 
is that the value of to where Re K t = Im K t is also the value 
of to where dlm/c t /dw = (i.e. where the phase is n/4, see 
Figure [3] final panel). It will turn out that this property is not 
always obeyed. 

4.2. Effects beyond FOSA and r approximation 

While a S function perturbation is disadvantageous numer- 
ically, it can be illuminating. If we impose on a flow with 



U = a test-field C with a 6(t) time dependence, then the 
value of c(0) depends only on the V • (uC) term in Equa- 
tion d2"Tl >. For f > 0, Equation (l2Tb reduces then to: 

^ = -V • (uc - uc) + kV 2 c (f > 0). (39) 
of 

Such a 6 perturbation then launches fluctuating fields which 
evolve according to an equation similar to Equation ( |39l . In 
passive scalar or kinematic dynamo cases, the evolution of 
the fluctuating field depends only on u, which is independent 
of the fluctuating field. The fluctuating fields will decay ex- 
ponentially according to turbulent or micro-physical diffusion 

or resistivity, but they will generate T (or £> in the magnetic 
case) for as long as they survive. It is the finite lifetime of the 
fluctuating fields that is at the physical core of this memory 

effect. 

In the passive scalar case we consider T = uc, to which 
the term V ■ (uc) in Equation (1391 does not contribute because 

mV ■ (uc) = 0. If the spatial dependence of our test scalar 
C is sinusoidal and lies only along a direction z, then c(0) 
will also have only sinusoidal behavior in that direction. If 
we imagine that the initial c(0) is proportional to sin kz, then, 
in the absence of other effects, two counter-propagating ver- 
tical streams with u z = ±u (assuming U — 0) will generate 
an advective sinusoidal signal from the V,(m~c) term of Equa- 
tion (f39): 



u z V z (u z c)(t, z) = 2m k cos kz cos coot, 



(40) 



where too = ku. In a turbulent system the above analysis can 
only be done for times shorter than or comparable to a turbu- 
lent correlation time t ~ 1/ku. For times larger than a turbu- 
lent correlation time, the standard e~^ T diffusion term will be 
important. A "turbulent" diffusion is formally possible even 
in steady flows, but it will be just the microscopic diffusion. 

We can combine the short timescale advective (oscillatory) 
and longer timescale diffusive (exponential) effects by a sim- 
ple multiplication: we expect the form for K t (t) to be similar 
to 

(41) 



K t (t) s u 2 0(f) e ,/T cos u> t 



Note that in a turbulent system we expect o>q ~ 1/t on di- 
mensional grounds and so the above analysis is not rigorous. 
However, as we will see in Section|6]this form fits the results 
reasonably well. In the Fourier space this becomes 



h(oS) 



1 



,2^2' 



(1 - iwr) 2 + cDqT 



(42) 



where £,o = tm 2 has been assumed, although this prefactor 
may not be accurate for ojq + 0. The corresponding Laplace 
transform is 

US) 1 f ST 



2^2- 



(1 + ST) 2 + COqT 



In the l imit ojq — > these expressions coincide with those of 
Section l4~Tl In Figure[3]we plot various representations of the 
integral kernel for different values of loot. 

In order to assess whether the proposed extension to cap- 
turing memory effects is viable, we shall use Equations d42l 
and d43l as fit formulae to determine the value of too and to 
find out how it depends on other aspects of the model such as 
the Peclet number and wavenumber of the mean concentra- 
tion. 
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In the absence of a detailed analogous motivation for a and 
77t we shall use in this paper Equations ( l42l i and (|43l l as fit 
formulae also in the magnetic case. In this case, we use these 
formulae for a and rj t and add corresponding subscripts a and 
rj to t and uq, where it replaces the subscript 0, i.e. we write 



a((jj) 



= A n 



1 



(1 - icOT a ) Z + (jJ l a T l a 



(44) 



Again, we expect co a T a and a> n T n to be of order unity, but in 
this paper we allow them to be adjustable parameters. Further, 
we use A a and A n as further fit parameters, modifying the 
amplitude. The relaxation times r and u> 1 and values derived 
from them such as a,o are merely characteristic times, and we 
do not attempt to laboriously average over the true values. 

Note that the above form for the kernel, Equation d4Tl ). is 
the simplest extension of the r approximation that qualita- 
tively fits our simulation results. From that perspective, we 



j/t(a>) 
fjtO 



1 - iu>T„ 



'(1 - iWT n ) 2 + Loftf 



(45) 
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Fig. 4. — Normalized time lag Ar/r versus lot for different values of loqt. 
Note the development of a peak near lot = 1 as loot is increased. 



replace Equation (l35l l by: 



2_2 



dt 



at 2 



at 



(46) 



Note also that, unlike Equation ( |36l >, for Equation ( l42l >. the 
value of w where the slope of the imaginary component is 
zero is not the same as the value of to where the phase is n/4 
(see the end of Section HTTt. 

As shown in Equation (IA5t of Appendix lAl for monochro- 
matic mean fields a phase shift <p K leads to a time lag 



At — <p K (to)/to, 



(47) 



so the flux T \j(t) depends only on the mean concentration gra- 
dient at time t - At and is given by -\k\ G(t - At). For the 
response function given by Equation (|42| >. the time lag is 



At 

T 



4> k {(jS) 



COT 



1 



-arctan 



1 + (to 1 - lo 2 } )t 2 

1 + (to 2 + tol)T 2 



(48) 



which always vanishes for large values of to and can have a 
peak near tor = 1 for too > with <x>*t =s 0.3273; see Fig- 
ured 



5. SIMULATIONS 

We consider two types of flows. For test purposes and com- 
parison with earlier work described in Section [2] we use the 
Roberts flow. The Roberts flow is given by 



with 



U = kflffZ - Z X Vl/r, 



ijt — (uo/ko) cos fcoxcos koy 



(49) 



(50) 



and kf = V2£o so that k 1 = k\ + ky, where k x — k y = ko is the 
wavenumber of the flow in the xy plane. This flow is capable 
of dynamo action once the magnetic Reynolds number, 



Re 



M 



(51) 



exceeds a critical value, Re^ > Re^crit = 5.52. [We recall 
that our test fields have spatial dependence given by k\, i.e. 
the smallest wavenumber that fits in the box. We note further 
that for ReM < 70 the most unstable wavenumber of the field 
that fits into the box is still \k z \ - k\, where k z was defined in 
Equation © and this agrees with the wavenumber of the test 
fields. Note however that, for Re^ = 100, for example, the 
most unstable mode would have \k z \ = 2k\.] 



The other alternative is forced turbulence. In that case we 
consider an isothermal equation of state with constant speed 
of sound, c s , and solve the momentum and continuity equa- 
tions 



dU 

at 



— = -U Vt7-c;fVlnp + /+p- 1 V-2pvS, (52) 



dp 

f t = -V • (Up), (53) 

where / is a random forcing function consisting of circularly 
polarized plane waves with positive helicity and random di- 
rection and phase, S is the traceless rate-of-strain tensor. The 
length of the wavevector of the forcing function, \kf\, is cho- 
sen to be in a narrow band around an average wavenumber 
£f. We adjust the strength of the forcing such that the flow 
remains clearly subsonic (mean Mach number is around 0.1). 
The details of the forcing function used in the present work 
can be found in Appendix A of Brandenburg & Subramanian 
(2005). For forced turbulence we define Re^ = u lms /T]kf. 

We consider a domain of size L x x L v x L z . In all cases, 
we take L x - L y - L z = 2n/ki. The ratio kf/ki is referred 
to as the scale separation ratio. Our model is characterized by 
the choice of fluid and magnetic Reynolds numbers as well 
as the Peclet number, based here on the wavenumber kf. The 
magnetic Reynolds number was defined in Equation (15 It , The 
fluid Reynolds and Peclet numbers are defined analogously, 



Re = u rm Jvkf, Pe = u Im jKk f 



(54) 



where the magnetic diffusivity 77 is replaced by the viscosity v 
and the molecular diffusivity k, respectively. 

We present the results in non-dimensional form by normal- 
izing k t (co), analogously to earlier work (Brandenburg et al. 
2008a), by 



l -TU 2 

2 rms 



(for the Roberts flow). (55) 



For turbulent flows, r is proportional to the turnover time, 
(Wrms&f) ■ However, in the limit of low Peclet number, mi- 
croscopic diffusion becomes important and dominates over 
the triple correlation terms. This means that the effective t 
is given by the microscopic diffusion time r = (Kk 2 y l . 

We define the Strouhal number as St = ru lms kf and can then 
write t as 

r = St/(u ms kf)- (56) 

The value of St characterizes the flow field. For turbulent 
flows of the form discussed in the present paper its value is 
of order unity (Brandenburg et al. 2004). Later in this paper 
we shall allow St to be a fit parameter. We present the results 
for a and rji by normalizing, depending on the nature of the 
flow with 



ao 



and 



a 



2 H rms ; 



3 "rms > 



t]t0 



-,-1 



2 "rms^f 



Mi-ms^f 1 (3D turbulence). (58) 



(Roberts flow) (57) 



3 i*rm\' x 1 



as was done in Brandenburg et al. (2008a). 

6. RESULTS 

Our choice of Equation (fTTT) results in transport coefficients 
that depend on the wavenumber of the mean fields. Through- 
out this section we will assume that our mean fields vary spa- 
tially according to k z = k\ unless otherwise specified. For 
simplicity therefore, we drop the fixed argument k z in T<{k z , to), 
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Fig. 5. — Real (solid circles) and imaginary (diamonds) parts of H(a>) for 
forced turbulence with kf/ki = 2.2, Re = 8, and Pe = 40. The solid and 
dash ed l ines are a fit using Equation {36} (with r determined using Equa- 
tion {56}, St = 2.7 as described in text). 



Fig. 6. — Real and imaginary parts of i(tS) for Pe = 100 for the Roberts 
flow. Note the perfect agreement with the fit formula {36} using r = 1/(kIc^) 
and K t o = ru 1 2 ms /2 (curves). 



a(k z , to) and ijt(k z , a>), and similarly for Tc{k z , s), a(k z , s), and 
fj t (k z , s). 

6.1. Passive scalar diffusion 

We now consider solutions of Equation ( T2lT i in the case of a 
turbulent flow, and consider first the case with a uniform gra- 

— — c0 

dient of C. This means that G is now constant in space, with 

G = Go(t). The resulting data agree well with the expres- 
sion Equation (l36b . where r is given by Equation d56l l with 
St = 2.7; see Figure[5] The fact that St > 1 should not be too 
surprising, because such a result has been obtained earlier for 
this flow, where r was estimated as the relaxation time in the 
t approximation (Brandenburg et al. 2004). 

The case of the Roberts flow, where U is obtained from 
Equations d49l ) and ( T50b is in some ways more interesting. In 

— c0 

the case of the same uniform gradient concentration C , the 
flux can be calculated analytically, as is done in Appendix |B1 
As the flow U is steady its correlation timescale is infinite 
and the only relevant relaxation timescale is the microscopic 
diffusion time t = (Kk^y 1 . The calculations result in the ex- 
pression K t (aj) = Kto/(l _ iwr); see Equation 1361 . This agrees 
with simulations as shown i n Fig ure [6] 

We suggested in Section 1431 that advective effects play a 
role only when the mean concentration gradient shows a vari- 
ation in some direction (i.e., a finite wavenumber), and we 
should not be surprised that Equation d36l l is adequate to ex- 
plain the transport of a passive scalar with zero wavenumber. 
The results of Figure [7] where a turbulent flow is used with 
Re = 8 and Pe = 40, and a sinusoidal variation of the mean 
concentration is imposed, are slightly better fitted with Equa- 
tion l42l than with Equation l36l l. 

The case of the flow U = uqzcosx with the same sinu- 
soidally varying concentration is discussed in Appendix ICl 
The value of to where Re/c, = ImAi> is not the same as the 
value of co where Im/c, has zero slope. Th is is implied by 
Equation (l36l >. as is discussed in Section |4~T1 In AppendixICl 
we present a simple one-dimensional model where the behav- 
ior is at odds with Equation (|36l l, although it can still be fitted 
reasonably well with Equation ( l42l . This example also illus- 
trates the difficulty in developing good and simple fits, as the 
fit parameters are expected to depend on the spatial variability 
(e.g. through too ~ ku). 
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lines are fits. 

6.2. Magnetic fields 

For small magnetic Reynolds numbers the functional forms 
of both a{a>) and ^ t (w) are similar to those in the passive scalar 
case. This is demonstrated here for the Roberts flow; see Fig- 
ure [8] where Re^ = 1, which is too small for dynamo action. 
Figure|9]shows the s dependence for the same case. However, 
for Re« = 10, which is large enough for dynamo action, the 
forms of a(u>) and ijtitj) look rather different; see Figure [TOl 
which is also for the Roberts flow. Qualitatively, the data are 
now closer to Equation (l42l . but a fit would be relatively poor. 
Therefore, we cannot rely on a fit to compute the correspond- 
ing Laplace-transformed kernel functions, which are shown in 
Figure[TT]for Re^ = 10 and 50. Note that, unlike the case of 
Figure|9]for Re^ = 1, for Re^ = 10 and 50, the slope of a(s) 
is positive. This is also a feature found by Hori & Yoshida 
(2008); see their Figure 10 for Re^ = 4, which corresponds 
to Re^ = 8 in our definition of the Roberts flow. 
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Fig. 8. — Real and imaginary parts of a(u>) and fj t (oj) for the Roberts flow 
with ReM = 1. The solid an d d ashed lines correspond to fits of the form 
Equation {36} using Equation {56} with St = 0.4. 

Figure [TT] allows us now to assess the error done by apply- 
ing the dispersion relation Equation ^ with constant values 
of a and T] t to cases where A + 0. A correct procedure would 
be to use a(s) and fjt(s) for s — A. This means that we must 
calculate 

AW^&Wh-iri+Ms)]^ (59) 

for s = A. These points can be obtained from the intersec- 
tion of A(s) with the diagonal, A(s) = s. For Re M = 10 
and 50 these values are at A{s = A) w 0.07M lms A:f and 
* 0.11 

^rms^f » 

respectively. By contrast, A(s — 0) ~ 0.04« lms £f 
and « 0.07« rms £f for these two values of Re^, respectively. 
These values are now in full agreement with those of A gIovlt h 
seen in Figure[T] This suggests that the reason for the discrep- 
ancy between the two curves in this figure is indeed connected 
with memory effects. 

Let us now turn to the calculation of a(to) and fj t (a>) in the 
case of turbulence. In this work we use kflk\ = 3, which 
is slightly larger than the values used earlier in the case of a 
passive scalar. This value is just large enough to allow for 
mean field dynamo action at the minimal wavenumber k - k\ 
(see Brandenburg et al. 2008c). For fef/fe] = 2.2 the scale 
separation between the scale of the forcing and that of the 
domain would be insufficient to allow for large-scale dynamo 
action (Haugen et al. 2004, Figure 23). 

By comparing runs of two different magnetic Reynolds 
numbers, Figure [12] for Re^ = 22 and Figure [T3] for Re^ = 
90, we can get some idea whether the features seen here are 
artifacts of small values of Re^f, or whether they begin to be 
of more general significance. The plots for aiu) and fjiioj) 
look similar and share the same basic features at both val- 




-0.4 



-0.2 



0.0 

ST 



0.2 



0.4 



Fig. 9. — Similar to Figure[8] but for the Laplace-transformed kernel func- 
tions a{s) and f) t (s) for the Roberts flow with Rem = 1. The fits are propor- 
tional to r/(l + st) and correspond to the fits used in Figure[8] 



TABLE 1 

Comparison of fit coefficients for a(a)) and fj t (ui) for forced turbulence. 



Re A a 



St„ 



St, 



22 
90 



1.00 
1.00 



2.00 
1.40 



1.20 
1.00 



0.48 
1.80 



1.40 
1.70 



0.55 
0.78 



ues of ReM, suggesting that the resulting fits for the response 
functions might be robust. We note that in all cases the phase 
shows a gradual transition from to n/2 as to increases, but 
it does not become negative (not shown). The corresponding 
fit parameters are summarized in Table [T] All the six non- 
dimensional fit parameters should be of order unity, and we 
see that this is indeed the case. Given that these values have 
unknown errors connected with the ambiguity in determining 
good fits, it is not possible to draw any serious conclusions 
from the trends that could be read off the table. 

Similar to the case of the Roberts flow, the fits to the Fourier 
transformed quantities are not perfect. Therefore we cannot 
use the Fourier transform fits to determine the corresponding 
Laplace transforms. In Figure [14] we show the directly deter- 
mined Laplace transformed values and compare with the fit 
inferred from Figure Q~3] However, in order to make the fits 
agree reasonably well, we have modified the amplitude fac- 
tors to A a = 1.37 and A n - 2.07. Note that the agreement is 
reasonably good, except near s = 0, where the actual growth 
rate is lower than what is inferred from the fit. This is related 
to the fact that the actual value of a(s) near s — is less than 
what is predicted by the fit formula. This suggests that the 
assumption of similar fit formulae both for a and rj t may be 
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Fig. 10. — Real and imaginary parts of a(ti)) and fjt(oj) for the Roberts flow 
with Rett = 10. Note that the lines are not analytical fits. 



too simplistic. 

As for the Roberts flow, the actual growth rate of the mean- 
field dynamo is obtained from the intersection with the diag- 
onal, which is shown as a dotted line in Figure [14] By solv- 
ing the induction equation for this flow for Re« = 90 we 
find that the actual growth rate is 0.04i< rrns A:f , which is clearly 
above the point where A(s) intersects with the diagonal (see 
the open symbol). However, this is to be expected, because 
for Rsm = 90 there is strong small-scale dynamo action so the 
actual growth rate will always exceed that expected from the 
mean-field dynamo. Such a discrepancy was noticed recently 
in connection with a study of the dependence of large-scale 
dynamo action on the magnetic Prandtl number (Brandenburg 
2009). 

7. DISCUSSION 

7.1. Frequency and growth rate dependence 

An important application of the present results is the deter- 
mination of dynamo growth rates. The usual dispersion rela- 
tion for isotropic helical turbulence predicts the growth rate to 
be 



A — ak - (77 + r]^)k (for constant a, rj t ). 



(60) 



However, if the resulting magnetic field really were to grow 
like e At , the effective values of a and rj t would be modified 
and would no longer be constant. By applying Equations (144-b 
and d45l l for a range of values of A for which 1 + At > we 
find that a and j] t become 



a(A) = a A a 



1 + Ar a 



(1 + AT a ) 2 + w2 T 2 



(61) 
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Fig. 11. — Laplace-transformed quantities, a(s), fjt(s), and A(s) for the 
Roberts flow with ReM = 10 and 50. Note the different signs of the slope 
at the intersection with the diagonal (denoted by circles). 



and 



1 +At, 



(1 + At % ) 2 + uItI 



(62) 



respectively. In these equations the occurrence of the terms 
<j,Tj for i — a or r] t is qualitatively new compared with earlier 
expectations based on the t approximation; see Section 14.11 
Note that the relaxation times r, and oscillation frequency <y, 
from Equation (l4"TT i are in general different for a and rj t ; see 
Table Q] 

A more direct way of calculating a(A) and rj t (A) is by using 
exponentially growing or decaying test functions proportional 
to e st , provided that 1 + st > 0, which sets the maximal decay 
rate for which equations (l6Tb and d62b are meaningful. The 
existence of a maximal decay rate is interesting: in such a 
system the fluctuating fields survive long enough to preserve 
the mean field. Clearly then, solutions of Equation (l60l , 

A = a(A)k - [tj + ri t (A)]k 2 (63) 

are required for self-consistent systems (be they dynamos or 
decaying mean fields). 

7.2. Wavenumber dependence 

In the work of Brandenburg et al. (2008a), which led to this 
paper, similar methods were used to find the dependences of 
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Fig. 12. — Real and imaginary parts of 5(ti>) and ?)t(<x> ) fo r turbulence at 
Re,w = 22. The lines denote fits to Equations {44) and (45) with A ff = 1, 
St tt = 2.0, (i) a T a = 1.2, and A, = 0.48, St, = 1.4, a> rl T = 0.55, respectively. 

a and 77, on the wavenumber k of the mean magnetic field. In 
that paper, it was shown for the Roberts flow that under FOS A 
we have 



a(k) 



1 + (a a k/kf) 2 



'/to 



1 + (a,A# f ) 2 



(64) 



where a a — a n = 1. They found that this result is also a 
good approximation to turbulent flows, but then a a and a n 
were treated as fit parameters that are of order unity. While 
that work noted that memory effects should be expected, they 
were not treated. Equation (|64i i can be directly compared to 
Equation d6Tb with the growth rate A set to which recaptures 
the test-field method as used in Brandenburg et al. (2008a). 
This might suggest that cos^of is related to the advection 
term in Equation fiTl i. so one might expect that cjq ~ ku ms . 
For a>o - St,A:Mrms then, the formulae from Equation (|6TT > and 
Equation d64b match exactly, and by capturing the dependency 
of or and rj t on past times, we are perforce treating the problem 
as also non-local in space. One might therefore be tempted to 
suggest that the combined dependence on u> and k could be of 
the form 



and 



a(k, to) = aoA a 



fj t (k, u) = 7] tQ A„ 



1 - ia>T n 



(1 



■ iojT a ) 2 + (a a k/k{) 2 ' 

1 - iwT, A 
(1-iwT,,) 2 + (a n k/k f ) 2 ' 



(65) 



(66) 



However, although such a formula is indeed obeyed in the two 
special cases a> — (Brandenburg et al. 2008a) and k — k\ 
(present work), some preliminary work suggests that this 
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Fig. 13. — Real and imaginary parts of a(aj) and fi t (w ) fo r turbulence at 
Re M = 90. The lines denote fits to Equations {44} and {45) with A a = 1, 



St tt = 1.4, a) a T a = 1, and A 



1.8, St 



: 0.78, respectively. 



equation is not valid in general, and that a multiplicative rela- 
tion of the form a(k, at) = a(k)a(a>) and rj t (k, co) = Tj t (k)fj t {a>) 
might be more accurate. 

7.3. Linear time dependence 

After our paper appeared as preprint (arXiv:081 1.256 lvl), 
Hughes & Proctor (2009) pointed out an inconsistency in 
the turbulent magnetic diffusivity tensor when allowing mean 
fields with a linear time dependence. They attributed this to 
the occurrence of a new contribution to the magnetic diffusiv- 
ity. In the following, we explain that their result is a natural 
consequence of using Equations (PHI and (|45l l, as advocated 
in our paper. 

The time dependence of the mean field in the paper by 
Hughes & Proctor (2009) is given by 



Bit) = B () + C t, 



(67) 



with constants Bo and Co. If we convolve this mean field with 
the kernels a and f), corresponding to the t approximation (i.e. 
proportional to e~^ Ta and e~'^ T \ respectively), we find the £ 
to be 

£(/) = (a Q - 7] t0 k){B + C t) - (T a a - T n J] t0 k)C - (68) 

This formulation matches the form of Equation (25) of 
Hughes & Proctor (2009), where their F is given by -(T a ao - 
T ;; ?7to^)- Re-expressing Equation d68l in terms of B(t) and 
dB/dt, as well as their curls, proportional J(t) and dj/dt, we 
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Fig. 14. — Laplace-transformed quantities, a(s), fj t (s), and A(s) for forced 
turbulence at Re^ = 90. In the last panel, the diagonal X = s is shown 
as a dotted line. The growth rate obtained by solving the three-dimensional 
induction equation, which allows for small-scale dynamo action, is indicated 
by an open symbol. 

can write Equation d68l l in the form 



p r T^r dB r dJ 

£ = a B - 77toy"oJ + r a — - V n ^—, 

at at 



(69) 



where T a = -u§T a and F, = -t]^T n quantify additional con- 
tributions to the mean electromotive force. In the more gen- 
eral case where cj a and (±> n are different from zero, we have 



1 



/ ,2-2 



-a T a 



(1 +<t£) 



2^2 \2' 



r,, = -TJoTr, 



1 



/ i2-r2 



(70) 



We recall that the formulation in Equation d69l only applies 
to the special case of variations of the mean field that are lin- 
ear in time. More generally, we have 



£(-D" 



17=0 



d"B 
~di" 



7 ( "W 



d n J 

dt" t 



(71) 



where 



a (n) = I a(t) f dt, T] ( " } = I 77(f) f dt. (72) 
Jo Jo 

These moments are related to the derivatives of a(a>) and fi t (a>) 
at <x> — with 



a 



in) = ( _ i)n ^ 

daj" 



dcu" 



(73) 



where the subscripts indicate that the derivatives are to be 
evaluated at (o = 0. Note, in particular, that F a = -o- (1) = 
Im(d5'/dw)o. 

Hughes & Proctor (2009) have computed values of r a for 
Re^f between 1 and 100 using a particular form of the mod- 
ulated wave flow of Otani (1993), referred to as MW+ flow, 
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Fig. 1 5 . — Real and imaginary parts of a(a>) for k = using the Otani ( 1 993) 
MW+ flow with ReM = 1 (upper panel) and Reu = 100 (lower panel). The 
normalization is chosen to be ag = uq. The insets show the scaling of Imu 
near the origin with slopes 0.2 and 7.8 for the upper and lower panels, in 
agreement with the results of Hughes & Proctor (2009). 

which is given by Equation d49l with k{ = kg and 

2wo 1 9 1 \ 

\p(x, y, f) = — — ( cos a>f t cos Icqx - sin Wf f cos koy J , (74) 

where a>f = uoko has been chosen. 

In order to substantiate our interpretation of their results we 
have computed a(io) for their case with k — 0. The result for 
the Fourier-transformed kernel is shown in Figure Q3] for the 
Otani MW+ flow with Re M = 1 and 100. Compared with Fig- 
ure [3] there are additional features related to resonances with 
the frequency of the Otani flow. Such features cannot be 
explained with our simple fit formula. This means that higher 
order terms will become important in those cases where the 
variation of the mean magnetic field is more complicated than 
just a linear increase. 

The value of T a can readily be read off as the slope of the 
graph of Im a> near the origin. Our results agree with those 
of Hughes & Proctor (2009), as is shown in the insets of Fig- 
ure Q2] We note, however, there are no good reasons to as- 
sociate the T a term with a correction to turbulent diffusion 
alone. Instead, there are corrections both to a and to rj t once 
the mean magnetic field shows strong time dependence. 

In this connection it is important to emphasize that these 
complications are mainly a consequence of the particular time 
dependence inherent to the Otani flow and are not typical of 
turbulence, as seen before. For Re^ = 100 there is a distinct 
spike at a> = a>f, while for Re^ = 1 there is a smaller spike 
at u> = 2u>f; see Figure [15] We hypothesize that these spikes 
are associated with the periodicity of the Otani flow. Simi- 
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lar behavior is known to occur for the Galloway & Proctor 
(1992) flow (Courvoisier et al. 2006), and is connected with 
the infinite correlation time of a flow with sinusoidal time de- 
pendence (Radler & Brandenburg 2009). 

8. CONCLUSIONS 

Naive application of the values of a and rj t to time- 
dependent problems can lead to errors. This is because 
the turbulent transport coefficients are in general frequency- 
dependent, due to memory effects. So, for each frequency 
and for each growth or decay rate (corresponding to imagi- 
nary frequencies) the transport coefficients need to be deter- 
mined separately. The full frequency dependence can then be 
used to calculate response functions via Fourier transforma- 
tion. The result can then be used to determine the response to 
general time dependences, including, for example, oscillatory 
growth. 

The response function formalism shows that one needs to 
know the past time history of the mean fields to compute tur- 
bulent transport correctly. This is not new, but what is new is 
the fact that the departures from the instantaneous approxima- 
tion can be quite substantial for flows such as the Roberts flow. 
For isotropic turbulence, on the other hand, the effects tend to 
be less dramatic and simple fit formulae with an exponential 
decay and an oscillatory part can be reasonably accurate. 

The presence of an oscillatory part in the response function 
proportional to cos uot leads to a sign reversal of a and %. 
Hori & Yoshida (2008) associate this with the "over-twisting" 
in illustrations of Parker's Q loops. In their picture, rising flux 
tubes may twist by more than 90°. This interpretation clarifies 
the naive expectation that a may change the sign when the 
Coriolis force becomes important. In fact, as our work now 
shows, such an effect would only occur if the mean magnetic 
field varies like a 5 function in time or if it shows other rapid 
variations. Conversely, for mean fields varying slowly in time 
the net a would not change the sign, although some past times 
are weighted negatively. 

In the present work, we have only looked at one type of 
memory effect, where the typical timescales in the integral 
kernel are comparable to the dynamical timescales of the tur- 
bulence. There is yet another type of memory effect that can 
occur on a resistive timescale, namely the one associated with 
magnetic helicity conservation. As explained in the appendix 



of Blackman & Brandenburg (2002), this is a purely nonlin- 
ear effect such that the relevant time scale becomes very long 
only when the magnetic field is strong. Obviously, this effect 
is not captured by our kinematic approach, nor would it be 
relevant in connection with the calculation of growth rates of 
the dynamo. 

The approach presented here may be useful for calculating 
memory effects of turbulent transport coefficients over a range 
of other related problems. Particularly important might be the 
question of the damping of acoustic waves by turbulent vis- 
cosity in the Sun, for example (Stix et al. 1993). Such damp- 
ing would lead to line broadening of the acoustic frequencies. 
The present work has demonstrated that such quantities can 
only be useful if one has a good idea of its frequency depen- 
dence relative to the frequency at which the turbulent viscosity 
is determined and the frequency at which it is to be applied. 

Our approach could also be useful in cases where the tur- 
bulence itself is time dependent. This would be relevant for 
modeling convection in pulsating stars. Such systems are cur- 
rently being treated with time-dependent mixing length the- 
ory (Gough 1977). It would seem appropriate to adopt inte- 
gral kernels also in that case. However, now there would be 
two frequencies to be considered: the frequency at which the 
turbulence varies and the frequency at which the mean field 
varies. Another problem is that the test-field method has only 
been used and tested in connection with magnetic and passive 
scalar diffusion problems, and has not yet been developed for 
calculating the components of the turbulent viscosity tensor. 
This would indeed be one of the outstanding problems in this 
field. 
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APPENDIX 

CONVOLUTION FOR MONOCHROMATIC VARIATIONS 



The purpose of this app endi x is to show that for monochromatic signals a convolution corresponds to a multiplication in real 
space. Consider Equation ( I A It for a monochromatic function 



G(f) = Gaj(t) = Go cos tot, 
where u> is a constant. Inserting this into Equation ( TTZb yields 



rut) 



-f 



k t (t - f)Go cos ojt f dt f = -G, 



k t (t - f )e dt' = -G Re e 



(Al) 



k t (t-t')e i0J{ '- n dt'. (A2) 



By using a change of variables one sees that the integral is just the Fourier transform of k t (t). Thus, we arrive at 

TJt) = -G Re [^^(w)] . (A3) 

The real part of k t shows therefore a modulation with cos a>t and the imaginary part with sin cot. By projecting against these two 
functions separately, we can determine the real and imaginary parts of k(co). Thus, the complex function k(co) can be obtained 

from TJS) as 

U.oS) = -2G0 1 {e [w, ¥ (1> (t)) t , (A4) 
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which is the result stated in Equation (f30b . The factor 2 stems from the fact that the average values of cos 2 cot and sin 2 cot are 1/2. 
This procedure can be trivially exte nded to tensorial relationships; cf. Equation d33l >. 
It is interesting to write Equation (IA3b by expressing K t {co) in terms of its modulus and its phase, |£ t | exp \<p K , so we have 



TJS) = -\k\ G(t - At), where At = tp K (co)/co, 



(A5) 



showing that memory effects change not only the amplitude of the effective transport coefficient, but they also lead to a time lag 
such that, for a given frequency, the mean flux is proportional to the mean fields at a certain later time. 

ROBERTS FLOW WITH OSCILLATORY MEAN CONCENTRATION GRADIENT 

As was generously pointed out by the referee, in the special case of a Roberts flow, Eqs. d49l ) and ( l50b , with a mean concentration 
C = zGq cos cot, we can solve the problem analytically. In this case, Equation ( |2"TT i becomes 



dc 2 
— = -it-Go cos cot - V ■ (uc - uc) + k\ c. 
at 

In a first step we employ FOSA and neglect V • (uc - uc), so the above reduces to 

kV 2 c = — V2moGo cos kax cos hoy cos cot. 

at 



This has as a solution 



where 



c(x, y,t) — - 



V2m ( )Go 



(co 2 + (Kkj) 2 )) 



— cos kox cos koy cos(cot 



Kk\ 



(c0 2 +(Kk^) 2 ) 1/r (CO 2 + (Kk 2 ) 2 ) 1 

Note that V • (uc - uc) = and so this particular solution is also valid beyond FOSA. We obtain then 



sin ( 



1 1/2 



T"(t) = — 



1S.G0 



1/2 



2(co 2 + (Kk 2 ) 2 ) 

We can now find the Fourier-transformed kernel through Equation (IA41 > 



cos(cot - cp)z. 



(Bl) 



(B2) 



(B3) 



(B4) 



(B5) 



k t (co) = -2G ' {e iM r (a (i)) t 



2{co 2 + (Kk 2 ) 2 ) 



1/2 



Kk 2 + ico 



(c0 2 + (Kk 2 ) 2 ) m ) 



\ 2 


1 + ICOTk 




' 1 ) 


) 


[l +co 2 t\) 


~ ~2~ 


^ 1 - \LOT k 1 



(B6) 



where we hav e defined r k 1 = Kk 2 . Equations d49l > and d50t imply that u\ = Mq/2, and Section l6Tl argues that t = t^. Accordingly, 
Equation ( IB6b reduces to Equation d36l l. 

SIMPLIFIED ONE-DIMENSIONAL MODEL 

A simple system that defies result (l36l l of the r approximation is one with a passive scalar whose concentration varies sinu- 
soidally along z with k z + and a steady flow u — (0, 0, u), such that u = u(x) = uq cos kox, so V • u — 0. The equation for the 
small-scale concentration then is 

dc — 1 

— = -V ■ (uC + uc-uc)+kV 2 c, (CI) 
of 



which becomes 



and in turn 



dc — dc dc 2 
— = —uG — u— + it — — + kV c, 

dt dz dz 



dc dc dc 2 n— 

— = -u— + u— + kW c - u(x)G(t,z). 

dt dz dz 

This system is linear, inhomogeneous, with variable coefficients. We note that G(t,z) - G(t)e' k ~ z , impose 

G(t) = Go cos cot, 

assume that 

c = c(t, x)e lkzZ + c.c. 



(C2) 

(C3) 

(C4) 
(C5) 
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and treat as the system as a time-dependent problem with complex c(t, x): 



dc ( '~ ' 

— = — \k z u{x) + kIc^ c + k — - + \k,uc — uqGq cos Icqx cos a>t . (C6) 
ot 1 <J ox 1 

Equation dC6b is the equivalent equation to Equation (IB2b (and reduces to that equation when k z = 0). The Fourier-transformed 
kernel can be calculated similar to Appendix [B] and in Figure [16] we present a numerical solution for the Fourier-transformed 
kernel for kjk\ = 5 and uo/iik, = 5. 
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Fig. 1 6. — Real (solid circles) and imaginary (open circles) components of K t for k. = 5 and u = ugz cos x (see Appendix|Cj> using as fit parameters A K = 0. 105, 
T a = 0.33, and a) K = 1.8. 
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